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Abstract 

We analyze special class of bipartite states - so called Bell diagonal states. In particular we provide 
new examples of bound entangled Bell diagonal states and construct the class of entanglement witnesses 
diagonal in the magic basis. 

1 Introduction 

In recent years, due to the rapid development of quantum information theory [1] the necessity of classifying 
entangled states as a physical resource is of primary importance. It is well known that it is extremely hard 
to check whether a given density matrix describing a quantum state of the composite system is separable 
or entangled. There are several operational criteria which enable one to detect quantum entanglement 
(see e.g. [2] for the recent review). The most famous Peres-Horodecki criterion is based on the partial 
transposition: if a state p is separable then its partial transposition (1(g) r)p is positive. States which 
are positive under partial transposition are called PPT states. Clearly each separable state is necessarily 
PPT but the converse is not true. We stress that it is easy to test wether a given state is PPT, however, 
there is no general methods to construct PPT states. 

In [3j (see also [I]) we proposed a class of bipartite states which is based on certain decomposition 
of the total Hilbert space C d ®C d into direct sum of ci-dimensional subspaces. This decomposition is 
controlled by some cyclic property, that is, knowing one subspace, say So, the remaining subspaces 
Si,... , E^-i are uniquely determined by applying a cyclic shift to elements from So- Now, we call a 
density matrix p a circulant state if p is a convex combination of density matrices supported on S Q . The 
crucial observation is that a partial transposition of the circulant state has again a circular structure 
corresponding to another direct sum decomposition So © ... ® S^-i- Interestingly, also realignment 
[5] leaves the circulant structure invariant. This class was generalized to multipartite systems [6]. Its 
separability properties were analyzed in [7j. 

The class of circulant states contains a subclass of states which are diagonal in the basis of generalized 
Bell states in C d <g>C d . The corresponding rank-1 projectors define d 2 -dimensional simplex known in the 
literature as magic simplex. Several properties of Bell diagonal states were analyzed [8j EJ HUJ [JTJ H2]- 
In the present paper we perform further studies of this special class of bipartite states. In particular we 
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provide new examples of bound entangled Bell diagonal states and analyzed the class of entanglement 
witnesses diagonal in the magic basis. 



2 Circulant states for two qudits 

Consider a class of states living in C d ® C d constructed as follows: let {eo, • • • , &d-i\ denotes an orthonor- 
mal basis in C d and let S : C d — > C d be a shift operator defined as follows 

Se k = e k+ i , (mod d) . (1) 

One introduces d d-dimensional subspaces in C d ® C d : 

S = span{e ® e Q , . . . , e d _i ® e d _i} , (2) 

and 

E n = (I®S n )S , n=l,...,d-l . (3) 

It is clear that S m and S n are mutually orthogonal for m ^ n and hence the collection {So, . . . , S^-i} 
defines direct sum decomposition of C d <g> C d 

c d ®c d = s e...es d _i . (4) 

To construct a circulant state corresponding to this decomposition let us introduce d positive d x d 
matrices = 0, 1, . . . , d — 1. Now, define <i positive operators p n supported on S ra via the 

following formula 

d-1 d-1 

Pn=J2 a if e v ® 5 " e « S * n = Yl a if ® e «+«J+« • ( 5 ) 

j,j=0 i,j=0 

Finally, we define the circulant density operator 

P = Po + Pi + • • • + Pd-i • (6) 
Normalization of p, that is, Tr p = 1, is equivalent to the following condition for matrices 

Tr (>) +a (l) + ... +a (<2-l) > ) = i . 



The crucial property of circulant states is based on the following observation [3]: the partially transposed 
circulant state p displays similar circulant structure, that is, 

(id <g> r)p = po e . . . p d _ x , (7) 

where the operators p n are supported on the new collection of subspaces S n which are defined as follows: 

S = span{e ® e^ 0) , e\ e 7r(1) , . . . , e d _i <g> e 7r(d _ 1) } , (8) 
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where n is a permutation defined by 

7r(jfe) = —k , (mod d) . (9) 

It means that 

tt(0) = 0, tt(1) =d - 1, ...,7r(d- 1) = 1 . 
The remaining subspaces £ n are defined by a cyclic shift 

S„ = (I®5 n )S , n=l,...,d-l. (10) 

Again, the collection {So, . . . , £<z-i} defines direct sum decomposition of C d ®C d 

C^C^So©...©^! . (11) 

Moreover, operators p n are defined as follows 

d-l d-l 

Pn = ®if ^ ® ^ e *(0w(j) ^ = "-if e *i ® e ^(i)+n l7 r(i)+n , (12) 

ij'=0 ij=0 

with 

= a (™+ m ) o (nS m ) , (mod d) , (13) 

m=0 

where II is a permutation matrix corresponding to tt, that is 

Hfc« = 5fc l7r (0 , (14) 

and o B denotes the Hadamard product of d x d matrices A and B. 

Theorem 1 A circulant state p is PPT iff > , for n = 0, 1, . . . , d — 1. 

It is clear that any circulant state p gives rise to the completely positive map A : M^(C) — >■ M^(C) 
defined as follows 

p=(id®A)P+, (15) 
where P^ denotes the maximally entangled state in C d ®C d , that is, 

^ d-l 

Pd = e ki®e kl . (16) 

k,l=0 

One easily finds the following formula for the action of A: 

d-l 



» 

n=0 



He kl ) = a kl e k+n,l+n • (17) 
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We call A a circulant quantum channel if A is unital, i.e. A(I) = I. It implies the following condition 
upon the collection of positive matrices a( n >: 

d-l d-l 

^2 ^2 O-tk e k+n,k+n = I • (18) 
fc=0n=0 

Note, that a dual map A# defined by Tr(pA(X)) = Tr(XA#(p)), is defined as follows 

d-l 

A # (ch) = ^2 a lk e k+n,l+n , (19) 
n=0 

i.e. it is defined by the collection of a^ T . It is well known that if A* is unital, then the original map A 
is trace preserving. Note, that in general A is neither unital nor trace preserving. In the next section we 
shall consider a special class of circulant states which give rise to unital and trace preserving circulant 
quantum channels. 



Example 1 A circulant state of 2 qubits has the following form 



boo 



no 



\ oio 



boi 



(20) 



an / 



where for a more transparent presentation we introduced matrices a := a' ' > and b := aW > 0. Note, 
that a circulant state (I3ip is usually called X-state in quantum optics community [14J. One easily finds 
for the partial transposition 



boo 



V 



'10 



aio 



aoi \ 

boi 



hi 



(21) 



where the matrices a = [aij] and b = [bij] read as follows 

~_ / O00 ^01 \ £ _ 

V ho au J ' 

Hence, p defined in is PPT iff 
and hence 



an I 



boo a oi 
«io 



a > and b > 



aoo a n > I 1 2 , ^oo°n > koi| 2 



(22) 
(23) 
(24) 
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Example 2 A circulant state of 2 qutrits has the following form 




where a := > 0, b := > and c := > 0. One easily finds for the partial transposition 

\ 



aoo 


^00 


coo 


&01 


coi 


aoi 


C02 


ao2 


bo2 




ho 


cio 


fell 


Cll 




C12 




bu 












an 




Ol2 








C20 




C21 




C22 






O20 


b20 




&21 




«21 




a 2 2 


b22 



(25) 



(26) 



aoo 


coi 


^02 


cio 


fell 


012 


b20 


021 


C22 



coo 


hi 


O02 


ho 


an 


Cl2 


CL20 


C21 


^22 



where the matrices a = [ay], b = [hj] and c = [cy] read as follows 

feoo a oi C02 

«10 Cn ^12 
C20 ^21 022 

For more examples see [3]. Interestingly, circulant structure is preserved under realignment. 

Proposition 1 The realignment of the circulant bipartite operator 

d-l 

A = ^ a^j Cij ® Ci+nj+n , 
n,i,j=0 



reads 



where 



d-l 



n,i,j=0 



(27) 



(28) 

(29) 
(30) 
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Example 3 The realignment of p denned in (|3ip leads to 



/ «oo 



R(p) 



OlO 



'All 



V b 



11 



0QO \ 



no 



an / 



Hence in this case one has 



( ooo 


boo \ 


, r {1) = { 


' aio 


&io A 


1 hi 


an J 


V °01 


aoi / 



Example 4 The realignment of p defined in (|25p leads to the circulant structure with 



aoo °oo coo 
en an 6n 

»22 C22 a22 



aio oio cio 
C21 a 2 i 621 
bo2 C02 ao2 



and it my be easily generalized arbitrary dimension d. 



a-20 D20 C20 
coi aoi 601 
012 C12 ai 2 



(31) 



(32) 



, (33) 



3 Generalized Bell diagonal states 

Consider now a simplex of Bell diagonal states [U [101 [Til 02] defined by 

d-l 

P — ^ j PmnPmn ? 
m,n=0 

where p mn > 0, ^ m ,n Pmn = 1 and 

p mn = (i ® c/ mn ) p+ (i c4 n ) , 

with C/ mn being the collection of d 2 unitary matrices defined as follows 

tt „ _ \mk on _ \mk 

with 

A = e 2m/d . 



(34) 



(35) 

(36) 
(37) 



The matrices U mn define an orthonormal basis in the space M^(C) of complex dx d matrices. One easily 
shows 

Tr([/ mn [4) = d (38) 

Some authors [13] call U mn generalized spin matrices since for d = 2 they reproduce standard Pauli 
matrices: 

Uoo = I , Uqi = o\ , Uiq = ia 2 , U u = 03 • (39) 
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Let us observe that Bell diagonal states (|34p are circulant states in C d ® C d . Indeed, maximally entangled 
projectors P mn are supported on X n , that is, 

n n = p 0n + . . . + , (40) 

defines a projector onto S n , i.e. 

s n = n„(c d ®c d ) . (4i) 

One easily shows that the corresponding matrices are given by 

Q W = HD^H* , (42) 
where -ff is a unitary d x d matrix defined by 

H kl :=^=X kl , (43) 

and Z)( ra ) is a collection of diagonal matrices defined by 

D$ := pfcnfo . (44) 

One has 



^ = - d EP^ m{k - l) , (45) 



m=0 

and hence it defines a circulant matrix 

4? = /£3 , (46) 
where the vector is the inverse of the discrete Fourier transform of Pmn 

(n is fixed). 

Consider now partial transposition of Bell diagonal states. One has the following 
Theorem 2 If d is odd all matrices are unitary equivalent 

~(n) = S n~(0) s U ^ (4?) 

for n = 0, 1, . . . , d — 1. If d is even one has two groups of unitary equivalent matrices: 

a^ = S k a^S^ k , (48) 

and 

for k = 0,1,..., d/2- 1. 
Therefore 

Corollary 1 Bell diagonal state is PPT if 
• a(°) >0ford odd, 



7 



(°) > and > for d even. 



• a 

The corresponding completely positive map A : M n (C) — > M n (C) is denned by the following Kraus 
representation 

d-i 

A(X) = £ p m nU mn XUl m , (50) 

m,n=0 

where p mn > and J2mnPmn = 1- One has A(I) = J2 mn p mn 1 = 1, which proves that A is unital. Note, 
that the dual map 

d-i 

A#(X) = PmnUl n XU mn , (51) 

m,n=0 

is unital as well. Hence, A defines unital and trace preserving quantum channel (doubly stochastic 
completely positive map). 

4 Special cases 

In this section we analyze special classes of Bell diagonal states. 

4.1 d = 2 

For 2-qubit case one obtains the following density operator 

a (n) = (Xn y n \ (52) 

\ Vn %n J 

where 

Xn = ^(POn+Pln) , Vn = ^{POn ~ Pin) , (53) 

for n = 0, 1. The state is PPT if and only if 

xl > M 2 , x\ > \y \ 2 . (54) 

The above conditions imply well known result that 2-qubit Bell diagonal state is PPT (and hence sepa- 
rable) if and only if 

Pmn < ~ ■ (55) 

4.2 d=3 

For d = 3 the Bell diagonal state is defined by the collection of 3 matrices 

%n %n %n 

a (n) = [ z n x n z n I , n = 0,l,2 , (5(i) 

z n %n X n 
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where 



and 



1 



;{POn +Pln+P2r 



(POn + Api n + Xp 2n ) 



Now, the PPT condition reduces to the positivity of 



Xq Z 2 Z\ 
Z 2 Xi z 
Z\ Z Q X 2 



> o , 



which is equivalent to the following conditions 



and 



xqX\ >\z 2 \ 



x xix 2 + 2Rez ziz 2 > x 1 1 2 + ^iki| 2 + x 2 \z 2 \ 2 



(57) 
(58) 

(59) 

(60) 
(61) 



Hence, even for d = 3 the PPT condition is by no means simple. It might considerably simplify if we 
specify x n and z n . Assume for example that zq = 0. Then ([60]) - ([6T]) imply 

(62) 



xqX\x 2 > xi\zi\ 2 + x 2 \z 2 \ 2 



4.3 d = 4 

For d = 4 the Bell diagonal state is defined by the collection of 4 matrices 



' Xn z n y n z n N 

Zn X n Z n Vn 

Hn Z n X n Z n 

\ Z n Vn Z n X n J 



n = 0,1,2,3 , 



(63) 



where 



Xn = J (POn + Pin + P2n + P3n) , 
Vn = J [POn ~ Pin + P2n ~ P3n) , 
Zn = jCPOn + Win ~ P2n ~ WSn) 



(64) 



Bell diagonal state of two qutrits is PPT iff 



r(o) 



/ x z 4 y 2 zi \ 

za x 2 zi y 

y 2 zi x z A 

\ zi y z 4 x 2 ) 



> , and a (1) 



/ xi z y 3 z 2 \ 

zo x 3 z 2 yi 

y 3 z 2 xi z 

\ z 2 yi z x 3 / 



> 



(65) 
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4.4 Special form of p mn 

Consider now special examples of Bell diagonal states by specifying the structure of probability distribu- 
tion Pmn- Let 

Pmn = 5mkK n , (66) 

with 7Tq + . . . + Ttd-i = 1. It gives rise to 



d-l 
n=0 



(67) 



For example if d = 2 and k = one obtains 



7Tl 



7Tl 



\ vr 



71"! 



(68) 



7TQ / 



This state is separable if and only if ttq = ttx = 1/2. One easily generalizes this observation as follows 



Proposition 2 i?eZZ diagonal state (61) is separable if and only if 

1 

7T = . . . = TT d -i = - . 

d 

Another characteristic class corresponds to 

Pmn — QmPn i 

i.e. p mn represents the product distribution. One has 

d-i 



P= PkiPki = A) © ■ ■ ■ © Pd- 



i - 



k,l=0 



where 



d-l 

Pn — Pn ^ Qm-Prr 



m=0 



Note, that matrices are related as follows 



a (n) = p n a 



where the matrix a reads 



d-l 

kl ~ 1 qri 



m=0 



Proposition 3 Bell diagonal state (11) is separable if and only if 



Po 



Pd - 1= d- 



(69) 
(70) 

(71) 

(72) 

(73) 
(74) 

(75) 
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4.5 Generalized lattice states 

Consider now a family of unitary operators acting on N copies of C d 

U m n = U mim (£> . . . ® U mNnN , (76) 

where m = (mi,...,mjv) and n = (m, . . . , rt/\r). It is clear that [7 mTl defines a family of D 2 = d 2N 
unitary operators in in C D = C d ® N . Note that 

(77) 

Now, let IV'J) denote a maximally entangled state in C D ®C D defined by 

= ^7= Yl ek ® efc ' ( 78 ) 

where 

efc = e fel <g> . . . <g> e fcjv . (79) 
One defines a family of maximally entangled states by 

|VW> = (I®U mn )W+) . (80) 

These states are parameterized a point (m, n) in the iV-dimensional lattice L^jL consisting of D 2 points. 
Now, a generalized lattice state is defined by a collection of points from L^y. for any subset / C ^ny) 



yr ^ ^mti , (81) 



Pi 

(m,n)£/ 

where P mn = \ipmn) (VWi| an d |I| stands for the cardinality of /. Clearly, 1 < |J| < |^|^| = D 2 . Let us 
observe that the above construction generalized a class of lattice states presented in |15t [TBI 02] • Lattice 
states of Benatti et. al. correspond to d = 2. In this case U mn are defined in terms of Pauli matrices (see 
formula (1391)1. 



5 Bound entangled Bell diagonal states 
5.1 Two qutrits 

Consider the following family of Bell diagonal states 

Pe = N e (P 00 + eUx + e" 1 ^) , e > , (82) 
where the projectors 11^ are defined in (I40p and the normalization factor reads 

N e = r . (83) 
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It corresponds to (cf. formulae (f6i| and (|58j) ) 



N e N £ N £ _! 

= -y , a?i = -g - e , x 2 = — e , (84) 

and 

2o = -y , zi = Z2 = . (85) 

Hence, conditions (|6U|) - (|61|) are trivially satisfied showing that (|82|) defines a family of PPT states. Now, 
it is well known [21] that p e is separable if and only if e = 1. Hence, for e ^ 1 it defines a family of bound 
entangled state in C 3 (8>C 3 . The entanglement of p e can be detected by using a realignment criterion 
[5]. In the next section we show that it can be detected also by the Bell diagonal entanglement witness. 
Note, that asymptotically 

lim Pe = o n 2 > lim = o n i ' ( 86 ) 

that is, one obtains separable states defined by normalized separable projectors onto £2 arid Ex, respec- 
tively. 

5.2 Two qudits 

Consider a family of states in C d ® C d defined by [181 120] 

7 »,i=o 



where d x d matrices 



4 7 



with 



dj for i 7^ j , 

e o + a 7 en + Y^iZv e « + &7 e d-M-i for i = j = , (88) 
S^A^Sti-i for i = j ^ 1 



a 7 = ^(7 2 + d-l), 6 7 = i( 7 - 2 + d-l) , (89) 



and the normalization factor reads 

s 

It gives the following spectral decomposition 

d-2 



7V 7 = d 2 - 2 + 7 2 + 7 ~ 2 . (90) 



£=2 
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In particular for d = 3 one obtains the following matrix representation: 



/ 1 



V 1 



1 \ 



6 7 



— (3P o + « 7 ni + 6 7 n 2 ) , (92) 



1 / 



with o 7 = |(7 2 + 2) , 6 7 = | (7 2 + 2) and the normalization factor iV 7 = 7 + 7 2 + 7 z 



6 Bell diagonal entanglement witnesses 

Interestingly many well known entanglement witnesses displaying circulant structure are Bell diagonal. 
It is well known that any entanglement witness W can be represented as a difference W = W+ — W-, 
where both W+ and W— are semi-positive operators in B(C d <S>C d ). However, there is no general method 
to recognize that W defined by W+ — W_ is indeed an EW. An interesting class of such witnesses may be 
constructed using their spectral properties [22l [23] . Let tp a (a = 0, 1, . . . , d 2 — 1) be an orthonormal basis 
in C d (&C d and denote by P a the corresponding projector P a = \ip a )(ip a \. Now, take d 2 semi-positive 
numbers X a > such that X a is strictly positive for a > L, and define 

L d 2 -l 
a=0 a=L 

where L is an arbitrary integer < L < d 2 — 1. This construction guarantees that W+ is strictly positive 
and all zero modes and strictly negative eigenvalues of W are encoded into W-. Consider normalized 
vector tp £ C d ®C d and let si(ip) > ... > Sd{ip) denote its Schmidt coefficients. For any 1 < k < d one 
defines A:-norm of ^ by the following formula 

k 

u\\l = Y, s2 M) ■ ( 94 ) 

It is clear that 

ll^lli < NV'lb < ... < HVIU • (95) 

Note that gives the maximal Schmidt coefficient of ip, whereas due to the normalization, WiftWj = 

= 1. One proves |22j the following 

Theorem 3 Let ^2^=0 HV'alll < 1 • -tf ^ e following spectral conditions are satisfied 

K>Hk, ol = L, . . . , d 2 - 1 , (96) 
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to 



where 

t/ien TV is an k-EW. If moreover X^q=i H^alljL-i < 1 anc ^ 

Mfc+l > A a , a = L, 

then W being k-EW is not (k + 

Let us observe that if is maximally entangled then 



,d 2 - 1 



Consider, therefore, the family of Bell diagonal states ip mn . On has the following 
Corollary 2 If L < d and 

A a > to , oi = L, . . . , d 2 — 1 , 
with \x\ = ^2a=o , i/ien PV^ = W+ — W- 7 - defines Bell diagonal entanglement witness. 

Example 5 Consider well known entanglement witness in d = 2 represented by the flip operate 



/ 1 



F 



\ 



1 / 



Note that 



F — Poo + Pio + -Poi — Pu , 
which proves that F is Bell diagonal and possesses single negative eigenvalue. 

Example 6 A family of EWs in C 3 ® C 3 defined by 





( a ■ ■ 
■ b ■ 


■ -1 • 


..-IX 




■ c 






W[a,b,c] = 


-1 • • 


c 

a 


■ • -1 




[-1 ■ ■ 


■ -1 • 


• c 

• • a J 



with a,b,c > 0. Necessary and sufficient conditions for W[a, b, c] to be an EW are 
1. < a < 2 , 
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2. a + b + c>2 , 

3. if a< 1 , then 6c > (1 - a) 2 . 

A family W[a, 6, c] generalizes celebrated Choi indecomposable witness corresponding to a 
c = 0. One finds the following spectral representation 

W[a, 6, c] = (a - 2)P 00 + (a + l)(Pio + Pjo) + bU t + C n 2 , 

which shows that W[a, 6, c] is Bell diagonal with a single negative eigenvalue 'a — 2'. 

Example 7 Consider a family of E Ws defined by 



6=1 and 



(104) 



W: 



\,fj, 



1 






-1 




• -1 


• l + n 


A 


A 


1* 




A • 




A 


A 






A • 


-1 






1 




• -1 


A* 






• 1 + n 


A» 




H 


A 


A 




l + /i 


A • 


-1 






-1 




• 1 



(105) 



with A, n > 0. Note that Wo,o = W[l, 1, 0]. One obtains for the spectral decomposition 

^A, M = -3P o + 2n + III + 3/iPoi + 3AP 02 • 
Let 7 > 0. One shows [20] that for 



A < 



and 



2 + 7" 2 ' 

l-7 2 -A(2 + 7 - 2 ) 
2 + 7 2 



Wa )A4 defines an indecomposable EW due to TV(W^ M p 7 ) < 0. 

Example 8 Entanglement witness corresponding to the reduction map A(X) = ITrX — X 
One has 

1 1 d ~ X 

W=jI®I-P+ = -Y,Pki-Poo, 



(106) 
(107) 
(108) 

in M d (C). 
(109) 



k.l=0 



which shows that W is Bell diagonal with a single negative eigenvalue (1 — d)/d. 
Example 9 A family of E Ws in C d ® C d defined bv pi [20] 



d-l 
i,j=0 



(110) 
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where the d x d matrices xf- k are defined as 




(d—k— l)eu + Eti S- 



en S l for i = j 
for i ^ j 



(111) 



It is well known |18| that Wd,k defines an indecomposable EW for k = 1, 2, . . . , d — 2. For k = d — 1 it 
reproduces the witness corresponding to the reduction map. Note that for d = 3 and k = 1 it reproduces 
W[a, b, c] with a = b = 1 and c = 0. On easily finds the following spectral representation 



showing that Wdk is Bell diagonal and the single negative eigenvalue corresponds to the maximally 
entangled state Poo- 

7 Conclusions 

We analyzed a class of bipartite circulant states which are diagonal with respect to generalized Bell 
(magic) basis. Such states are characterized by an elegant symmetry which considerably simplifies their 
analysis. We analyzed several examples of bound entangled states and provided corresponding entangle- 
ment witnesses which are Bell diagonal. 
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